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We describe a generalization of the Lellouch-Luscher formula to the case of multiple strongly- 



Q_i coupled decay channels. As in the original formula, our final result is a relation between weak 

matrix elements in finite and infinite volumes. Our extension is limited to final states with two 



scalar particles, with center of mass energies below the lowest three- or four-particle threshold. 
Otherwise the extension is general, accommodating any number of channels, arbitrary strong 
coupling between channels, as well as any form of weak decay operators in the matrix elements. 
Among many possible applications, we emphasize that this is a necessary first step on the way to 
a lattice-QCD calculation of weak decay rates for D — > %% and D — > KK. Our results allow for 
arbitrary total momentum and hold for degenerate or non-degenerate particles. 
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1. Introduction 

Lattice calculations of weak decay amplitudes for processes with multiple hadrons in the final 
state are more difficult than calculations for leptonic or semi-leptonic decays. One of the many 
challenges faced is that finite volume effects are more important in the former case. 

For the simplest hadronic decay, K — > KK, the finite volume effects were worked out by 
Luscher and Lellouch over a decade ago [1-5]. This work, together with a great deal of addi- 
tional development in both formalism and computational techniques, has allowed for significant 
progress towards a first principles determination of these weak matrix elements [6, 7]. Specifically, 
results for the isospin-two sector with a complete error budget are now available [7], and the more 
challenging isospin-zero sector is expected to follow within the next few years. 

We are thus led to consider what information lattice calculations might eventually offer con- 
cerning heavier meson decays. For example, LHCb and CDF each recently reported CP-violation 
in the difference of CP-asymmetries for D° — > n~7t + and D° — > K K + [8, 9]. Lattice simulation is 
the only known method to check consistency of these results with the Standard Model. Therefore, 
though such a calculation is clearly very challenging, it is natural to begin investigating how it 
might be done. In particular, in contrast to the K nn case, the finite volume effects for decays of 
heavier mesons are not yet well understood. Our aim is to take a first step towards improving this 
understanding. Details can be found in the longer publication, Ref. [10]. 

Specifically we show that, if one can ignore all but two-particle channels, then a generalization 
of the work of Lellouch and Luscher allows for a first principles calculation of D decay amplitudes 
from the finite volume spectrum obtained via lattice simulation. Incorporating multiple strongly 
coupled two-particle channels is a precursor to the more difficult challenge of also coupling in four 
and higher particle states. Indeed such states, for example four pions, are expected to be important 
at the D mass (M D o = 1865MeV). 

Our discussion begins with a statement of the multiple-channel quantization condition. This 
was obtained using a field theoretic derivation in Ref. [10] and is a necessary first step before de- 
riving the Lellouch-Luscher formula. We note that there have been a large number of recent papers 
studying the generalization of the Luscher quantization condition to multiple two-body channels 
and assessing its utility [1 1-15]. After giving the quantization condition, we then show how it may 
be used to derive the generalized Lellouch-Luscher formula, which allows one to extract amplitudes 
for decays into multiple, strongly coupled, two-particle states. 

2. Generalized LUscher quantization condition 

In this section we state, without derivation, our extension of Liischer's quantization condition 
[1^1], to the case of multiple, strongly-coupled channels. As emphasized above, the extension 
allows for any number of coupled channels, each containing two scalar particles. 

Throughout the article we take finite volume to mean a finite, cubic spatial volume of extent 
L with periodic boundary conditions. We take L large enough so that exponentially suppressed 
corrections can be ignored, and take the (Minkowski) time direction to be infinite. 1 Although 

'Lattice simulations are performed with Euclidean time, and the spectrum is found from the exponential decay of 
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lattice simulation is the target of our investigation, we assume here that discretization errors are 
small and controlled and therefore use continuum field theory (zero lattice spacing). 

The kinematic variables which enter our discussion are the total energy E, the total momentum 

p = 2 _^l_ ^ p£ Z 3 ), (2.1) 

and the center of mass (CM) frame energy E* = \/ E 2 — P 2 . The condition that only two-particle 
channels are open is effected by requiring that E* sit below any higher particle state that is coupled 
to the sector of interest. In the case that a symmetry prevents even/odd coupling, this means taking 
E* below the lowest four-particle threshold. 

The quantization condition, which we derive in Ref. [10] by following Ref. [16], is of the form 

A-*{L,E*,P)=0. (2.2) 

Here j$ represents the infinite volume scattering amplitudes and is defined explicitly in Eq. (2.3) 
below. The explicit form of A"^ is also given later in this section [Eq. (2.7)]. To understand the 
utility of Eq. (2.2) one should first suppose a theory of scalar particles, for which all two-to-two 
scattering amplitudes are known. Then, at fixed {L,P}, becomes a known function of a single 
variable, E*. The content of Eq. (2.2) is that one should search this function for all roots, 
(k = 1 , 2, 3, • • •)■ These then give the energy spectrum of the finite volume theory. 

We comment that this description is actually the reverse of the most common case, in which 
one has energy levels from a lattice simulation and wants to determine physical scattering ampli- 
tudes. In this case one parametrizes the scattering amplitudes at fixed energy, Eq, with a number K 
of real unknowns. Next one must determine K sets {Li,Pi},- • -{Lk,Pk} which put £g in the sim- 
ulated spectrum. In this way one can produce K independent equations with the form of Eq. (2.2) 
and solve for the scattering amplitude. At the end of this section, we describe this method in detail 
for a particular toy model. 

We now turn to a precise description of the scattering amplitudes j$ . We assume a total of 
N open two-particle channels, labeled i,j = 1, • • • ,N. Then the two-to-two scattering amplitudes at 
fixed E* depend only on (a) the particle content of the in and out-state and (b) the relative direction 
of motion between incoming and outgoing particles in the CM frame. It is convenient to give the 
scattering amplitude redundant dependence by including both incoming and outgoing CM frame 
directions of motion, labeled k* and k*' respectively. The amplitude may then be written 



^ij(k*,k* ) = AltJKij. Am . l2jm2 Y hm {k*)Y^ nn {k* ) , (2.3) 

where a sum over i\,m\,l2,m2 is implied. We comment that ^y^ 1)OTl ^ 2 m2 is not diagonal in 
channel (if) space but is diagonal in angular momentum (£±, m\\ lt^m%) space: 

°^ij-A m ;t 2 ,m 2 = ,m §hh 8m im2 (2-4) 
(no sum). This is due to the rotational invariance of the infinite volume theory. 



correlators. In the corresponding Minkowski theory, the same spectrum is given by energy poles in momentum-space 
correlators. For our discussion the latter description turns out to be more convenient. 
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As a specific example, consider a toy model with three types of scalar particles: pions (with 
mass Mjc) and kaons and anti-kaons (degenerate with mass M£). Suppose, as is true for the isospin 
zero sector of physical pions and kaons, that the two pion (Tin) and kaon anti-kaon (KK) states 
only couple strongly to each other and also to states with four or more particles. It follows that, if 
we assume unphysical masses which satisfy 2Mk < AM n and we also require that E* lies below the 
four pion threshold, then the only open channels are kti and KK. In this case N = 2 with i = 1 for 
the pions and i = 2 for the kaons. Then, for example, ^#12 is the KK — > %% scattering amplitude. 

Returning to arbitrary two-particle channels, the infinite dimensional matrix ^ij-/ umi ^ 2 , m2 is 
one of the two ingredients needed to define A . The other is a kinematically determined matrix F 
which acts on the same space and is given by 



Rexj i 
-oe l e 2 o mim2 + 



ALE* 



(2.5) 



where x ; is the solution to 
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(2.6) 



with Mj i and My 2 equal to the two-particle masses of the jth channel. 3ff m is a generalized zeta 
function defined in Ref. [17] and T] = 1/2 for identical and 1 for non-identical particles. We com- 
ment that, in contrast to ^ , F is diagonal in channel space but not diagonal in angular momentum 
space. This is due to the breaking of rotational symmetry by the finite volume condition. 
We now state, without proof, the explicit form of our secular equation (2.2): 

A-*(L,E*,P) = det(F- 1 



+ 1 



0. 



(2.7) 



This result agrees with that found in Ref. [18]. It is also agrees with the earlier work of Ref. [13] 
in the limiting case of only s- wave scattering and P = 0. 

We next note that Eq. (2.7) is only practically useful in the case that the scattering amplitude is 
negligible above some £ max (^^ max,m =0). In this case one can show, by extension of a proof in 
Ref. [16], that F may also be truncated with no additional approximation. This is nontrivial since 
F is not diagonal in angular momentum. However, the projection in ^# turns out to be enough. 

We now return to our toy model of pions and kaons and incorporate the additional assump- 
tion that the s-wave scattering amplitudes dominate (£ max = 0). In this case Eq. (2.7) takes on a 
particularly simple form: 

^ • 



A^(L,E*,P) = det 







[F{(L,E*,P) 



'1^2 



E* 



^U(E*) 



= 0, 
(2.8) 



where all entries are now numbers and all functional dependence has been made explicit. 

We conclude this section by describing the extraction of scattering amplitudes in the s-wave 
TiTi-KK model. Due to unitarity and symmetry of the 5-matrix, the scattering amplitudes for this 
model are given at any fixed CM energy, Eq , by three real parameters. For an explicit parametriza- 
tion see Ref. [10]. One must therefore determine, from lattice simulation, three sets {L\ ,Pi}, {Z^,^}, {^3,-^3} 
which put E^ in the spectrum. These may then be substituted into Eq. (2.8) to deduce three inde- 
pendent equations, which constrain the three unknowns and thus determine ^-^^Eq). 
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iAMjk = 



Figure 1: The diagram giving rise to the amplitude perturbation t^M . 

3. Generalized Lellouch-Liischer formula 

In this section we show how the quantization condition of the previous section may be used to 
derive the generalization of the Lellouch-Liischer formula, relating weak matrix elements in finite 
and infinite volume. We restrict ourselves, for the entire section, to the s-wave %%-KK model intro- 
duced above, and therefore only reference the specific quantization condition (2.8). The arguments 
given here can be easily generalized to any number of two-scalar-particle states. 

We begin with a precise statement of the problem. We introduce an operator, J%w(x), which 
weakly couples the %% and KK to an otherwise non-interacting one particle state, which we sug- 
gestively label D. Our aim is to derive a formula which takes as input ^# and also the finite volume 
matrix elements 

M D ^ n =(n\je w (0)\D), (3.1) 
and gives from this the infinite volume decay amplitudes 

A D ^ z = (KK\jr w {0)\D) A D ^ K x = {KK\Mb{G)\D). (3.2) 

Because our quantization condition is only valid below the four-pion threshold, we require Md < 
4M n , where Mo is the D mass. This, together with 2Mk < 4M n , means that our result is not an 
honest description of physical D decay. What we present here is a necessary first step towards that 
challenge. 

The derivation proceeds by incorporating the weak Hamiltonian into the original theory via 

Jf?(x) — > J^(x)+XJ^ w (x), (3.3) 

where A is a real parameter that we can freely vary and, in particular, make arbitrarily small. The 
new Hamiltonian changes the scattering amplitudes and therefore also changes the form of the 
quantization condition: 

A-* +A ^{L,E*,P) = 0. (3.4) 
The leading order (in A) form of the change to the amplitude t^M is proportional to the matrix 

This change is due to the diagram shown in Figure 1 in which the two scalar particles of the in-state 
combine to form a virtual D which then decays into the two scalar particles of the out-state. 
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We now suppose that {L^P^} are chosen so that, in the absence of any weak interaction, Md 
is in the two-particle spectrum: 

A^(L k ,E* =M D ,P k )=0. (3.5) 

Here k = 1 , 2, • • • labels the different solutions to the quantization condition at E* = Md- Now comes 
the key point: When the modified Hamiltonian [Eq. (3.3)] is used, then the values L^,Mo,P^ no 
longer satisfy the quantization condition. However, by using leading order degenerate perturbation 
theory, one can correct the energies to restore the condition's validity. We conclude that 

A^ +A ^(L n ,E* =M D ±lL 3 \M D ^ n \,P n ) = 0. (3.6) 

Actually Eq. (3.6) is only valid through order A. Indeed, it is by expanding in X and requiring 
that the linear coefficient vanish that we reach our generalization of the Lellouch-Liischer formula. 
The final result may be written as 

\Cji(L n ,P n )AD^7C7c +CK[L n ,P n )A D ^ fK g\ = \Mo^ n \ ■ (3.7) 

This is the main result of this section. The coefficients C % and Ck depend on Md, {L n ,P n }, ^(E* = 
M D ) and dJt{E* = M D ) jdE*. We direct the reader to Ref. [10] for their specific forms. 2 

To understand the use of Eq. (3.7) one need only recall that the quantization condition Eq. (2.8) 
may be used to determine the scattering amplitude in the region of E* = Md- To get the amplitude 
one needs at least three sets {L\,P\}, {L2,P2}> {£3,^3} which put Md in the spectrum. Indeed 
one also needs three sets which put a slightly higher energy in the spectrum, in order to extract 
the derivative. Since three volume-momentum values are needed anyway, it is natural to calculate 
the corresponding three values of |M/)^„| (by lattice simulation) and also the coefficients C % and 
Ck (by substituting determined quantities into the equations of Ref. [10]). This results in three 
independent equations of the form 

I (known number )Ad^,kk + (known number )A D ^ K g\ = (known number) . (3.8) 

Indeed it turns out that three independent equations are precisely what is needed to determine 
the weak decay amplitudes. To see this one must first show, as we do in Ref. [10], that Ad^hk and 
Ad->kk ma y b e written as a known linear combination of two real numbers, which we call vi and 
V2? The three equations (3.8) may therefore be recast in the form 

I (known number)vi + (known number)v2| = (known number) . (3.9) 

Any two equations constrain the values of the real parameters up to sign ambiguities which are 
lifted by the third equation. 

2 The notation of [10] is slightly different, and is related to the notation in this article by 

C x = \tf M 2\- 1 / 2 Jqfr 1 [c 1 e- iS *c e -c 2 e- iS l>s E \ , 

Ck = \^r l/2 \ftlm [cie- i5a s £ +c 2 e- i5 P C£ \ . 

All parameters appearing on the right hand side are defined in [10]. Specifically "^ M i is defined in Eq. (88), q* in Eq. (34), 
c\ and C2 in Eq. (91) and S a , dp, c £ and s £ in Eq. (48). 

3 For weak operators that violate time reversal (T) symmetry there are some subtleties at this stage. We describe 
how to handle such operators in Ref. [10]. 



6 



Multiple-channel generalization of Lellouch-Liischer formula 



Maxwell T. Hansen 



4. Conclusion 

We have presented the generalization of the Liischer quantization condition to an arbitrary 
number of strongly coupled, two-scalar-particle channels in a moving frame. We have also sketched 
the generalization of the Lellouch-Liischer formula, which gives decay amplitudes into these chan- 
nels. Work is underway to further generalize this to final states with higher particle numbers. 
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